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ABSTRACT 

Global Positioning System (GPS)-based navigation has become common for low-Earth orbit 
spacecraft as the signal environment is similar to that on the Earth's surface. The situation changes 
abruptly, however, for spacecraft whose orbital altitudes exceed that of the GPS constellation. 
Visibility is dramatically reduced and signals that are present may be very weak and more 
susceptible to interference. GPS receivers effective at these altitudes require increased sensitivity, 
which often requires a high-gain antenna. Pointing such an antenna can pose a challenge. One 
efficient approach to mitigate these problems is the use of a digitally steered antenna array. Such 
an antenna can optimally allocate gain toward desired signal sources and away from interferers. 

This paper presents preliminary results in the development and test of a digitally steered 
antenna array for the Navigator GPS research program at NASA's Goddard Space Flight Center. 
In particular, this paper highlights the development of an array and front-end electronics, 
the development and test of a real-time software GPS receiver, and implementation of three 
beamforming methods for combining the signals from the array. Additionally, this paper discusses 
the development of a GPS signal simulator which produces digital samples of the GPS L1C/A 
signals as they would be received by an arbitrary antenna array configuration. The simulator 
models transmitter and receiver dynamics, near-far and multipath interference, and has been a 
critical component in both the development and test of the GPS receiver. 

The GPS receiver system was tested with real and simulated GPS signals. Preliminary results 
show that performance improvement was achieved in both the weak signal and interference 
environments, matching analytical predictions. This paper summarizes our initial findings and 
discusses the advantages and limitations of the antenna array and the various beamforming 
methods. 

INTRODUCTION 

The Global Positioning System (GPS) was designed for users on or near the surface of the Earth, 
and the GPS satellites correspondingly direct their transmissions toward the Earth. Space-based 
GPS receivers must deal with the high dynamics of orbital motion and the challenges of designing 
electronics suitable for the space environment, but GPS signal strength and availability in Low 
Earth Orbit (LEO) is similar to that on Earth. Indeed, GPS navigation for LEO spacecraft has 
become commonplace [1-5]. The situation is different when the receiver travels in high-altitude 
orbits above the GPS constellation where GPS signals become sparse and weak. Taking advantage 
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of signals from the sidelobes of the GPS transmitters can greatly improve signal availability above 
the constellation, but requires higher sensitivity in the receiver system [6-8]. 

One approach to increase sensitivity is through the use of special signal processing techniques 
for weak signals. Significant gains can be achieved in this way, but such weak-signal processing 
algorithms are often associated with reduced robustness against interference and degraded 
measurement quality. A complementary and more direct way to increase sensitivity is to gather 
more signal energy with a high-gain directional antenna. However, if the signals of interest arrive 
from a wide range of directions, as occurs in the case of GPS, physically pointing a high-gain 
antenna toward each desired direction becomes impractical. An equivalent result can be achieved 
with an electronically steered antenna array that does not require a mechanically gimbaled 
antenna. Such arrays offer advantages beyond simple single-lobe directivity in that they provide 
the potential to synthesize multiple beams in which gain is optimally distributed toward desired 
signals and away from interferers. 

This paper presents preliminary results in the development and test of a digitally-steered antenna 
array for the high-altitude Navigator space GPS receiver program at the Goddard Space Flight 
Center (GSFC) of the National Aeronautics and Space Administration (NASA). We refer to this 
project as the Navigator GPS Digital Beam Former (NavDBF). The next section gives a background 
about the Navigator program, followed by a brief review of antenna arrays and beamforming. 
Next, the GPS signal simulator and the three beamforming methods implemented in this project 
(a traditional deterministic beamformer and two adaptive optimization-based methods) are 
described in detail. The results of beamforming simulations and tests are presented and discussed. 
A summary of the project and possible future work concludes this paper. 

NASA GSFC'S NAVIGATOR GPS RECEIVER PROGRAM 

NASA GSFC's Components and Hardware Systems branch has been developing GPS receivers 
for space since the mid 1990s. Early receivers used commercial GPS chipsets and were tested on 
balloon experiments and in space-GPS research [6], including extending space-based GPS to high 
altitude missions above the GPS constellation itself [6, 8, 9]. At this altitude, signals are sparse 
and weak (see Figure 1). This research suggested the need for a specialized receiver designed with 
high-altitude applications in mind, and so was born the Navigator GPS program. 

The current-generation Navigator receiver, which achieves its high sensitivity through signal 
processing techniques, provides excellent performance in many high altitude applications [7, 8]. 
This capability led NASA to select Navigator as the onboard GPS receiver for the highly-elliptic 
Magnetospheric Multiscale (MMS) mission [8, 10], for which we have built and delivered eight 
flight-qualified receivers. Besides MMS, Navigator receivers have been built for the Global 
Precipitation Measurement (GPM) mission and an orbital demonstration during the Hubble 
Space Telescope Servicing Mission 4 (HST-SM4). The weak-signal capability of the Navigator GPS 
receiver has led to its use in numerous mission studies, such as the Geostationary Operational 
Environmental Satellite (GOES) and Orion reentry. The Navigator technology has also been 
commercialized and is used in the receiver for the Orion Crew Return Vehicle and a GPS receiver 
product developed by Broad Reach Engineering. 

NASA has continued to increase the performance of the Navigator GPS receiver through Internal 
Research and Development (IRAD) funding. The next-generation Navigator receiver will increase 
sensitivity by at least a factor of 10 to about 15 dB-Hz, improve performance in the presence 
of interference, and include support for modern GPS signals such as L2C and L5. The use of a 
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Figure 1: Regions of visibility of the GPS satellites. The dashed lines represent the approximate limits below 
which the side and main beams of a GPS satellite can be acquired by a GPS receiver sensitive to 25 dB-Hz 
(current Navigator capability) and 15 dB-Hz (proposed next-generation Navigator). The Geosynchronous 
Earth Orbit (GEO) and the two phases of the MMS mission orbit are shown for comparison. 


digitally-steered antenna array attempts to meet the first two of the these three goals and is the 
focus of the NavDBF project. 


Navigator Technology Development and Implementation Plan 

The Navigator GPS team has developed a three-stage technology advancement process to mature 
future Navigator GPS technologies from concept to flight (see Figure 2). This process defines a 
chain of three hardware platforms, which interface to low-level signal processing routines via 
an abstracted hardware layer and allows the flight software (with MMS and GPM Technology 
Readiness Level (TRL) 8 heritage) to be common among all three development stages. 


Stage 1. NavSDR (TRL 3/4) 

r \ 


Stage 2. NavDev (TRL 5) 

r \ 


Stage 3. SpaceCube (TRL 6/7) 
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Figure 2: Three-stage navigator GPS technology development process 

In the first stage of development, new algorithms and other technologies are implemented in 
our PC-based real-time Software Defined Radio (SDR) called the Navigator SDR (NavSDR). The 
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NavSDR provides flexibility and ease in developing and testing new algorithms, but is furthest 
from being flight-ready. Making the transition from stage one to two involves porting signal 
processing enhancements from the NavSDR software to the Very High Speed Integrated Circuits 
(VHSIC) Hardware Description Language (VHDL) for implementation on a Field Programmable 
Gate Array (FPGA). Our second stage platform, which we call the Navigator GPS Development 
Platform (NavDev) is based on a commercial FPGA development board (Xilinx ML605) and uses 
flight heritage firmware from the FPGA-based Navigator GPS receiver. This stage provides a 
low-cost, high-performance platform that is a short step away from flight. The third and final stage 
is a Xilinx Virtex 5 based flight avionics platform (GSFC's SpaceCube). The transition from stage 
two to three involves re-targeting the firmware for a slightly different FPGA family and replacing 
the laboratory Radio Frequency (RF) front-end hardware with flight electronics. The three stages 
of development correspond, roughly, to TRL levels 3-4, 5, and 6-7, respectively [11]. 

The plan for the development of the NavDBF was to build a 16-element antenna array 1 and 
front-end and implement digital beamforming algorithms in the stage one NavSDR platform. 
We planned to test the receiver with both live-sky and simulated GPS signals. A parallel effort 
to develop an array signal simulator in MATLAB, to produce simulated sampled data from an 
array and array front-end, was undertaken in order to test our beamforming algorithms before 
the fabrication of the front-end was completed. 

ANTENNA ARRAYS AND BEAMFORMING 

Antenna arrays have been extensively used in communication and radar systems that may be 
subject to intentional and unintentional interference. Adaptive antenna arrays were introduced in 
the early 1960s to obtain antenna directivity (see [12, 13]). The adaptivity in these antennas was 
achieved using a set of phase-locked loops to adjust the phases of each of the antenna elements 
relative to a reference signal. 

There are many different algorithms that have been used for adaptive electronic steering of 
antenna arrays. Two popular algorithms are the Least Mean Squares (LMS) and Minimum 
Variance Distortionless Response (MVDR) methods. The adaptive LMS method was introduced 
by Widrow and Hoff in 1960 [14] and subsequently applied to antenna arrays by Widrow in [15] 
and has been extensively used in various applications including radar, communication, and 
control systems. The MVDR adaptive algorithm was introduced by Capon in 1969 [16] where it 
was used to combine the output of an array of sensors for seismic applications. Both the LMS 
and MVDR adaptive algorithms have been, and continue to be, extensively used in various 
applications [17, 18]. In particular, for applications where the direction of the desired signal 
is not known a priori and the characteristics of the noise and the direction of interference are 
not completely known, adaptive methods have the potential to be far superior to conventional 
methods that produce fixed antenna patterns [17, 18]. 

Array signal processing has been applied to GPS receivers for interference mitigation (e.g., 
[19-28]). Moelker et al. [19], for example, considered adaptive antenna arrays and compared 
various algorithms and array structures. In [21], Brown and Gerein demonstrated sensitivity gain 
and interference mitigation in a research GPS array receiver. Lorenz and Boyd [22] demonstrated 
that a GPS antenna array can be used to mitigate multipath and interference errors using an 
algorithm based on angle of arrival differences. Kalyanaraman and Braasch [27] integrated an 

Hhe choice of a 16-element array provides the next-generation Navigator's target of 10 dB sensitivity gain, and 
more practically, was the largest array we could afford to build at the time. 
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adaptive array with a software GPS receiver and studied the effect of narrow band as well as 
wide band interferences on code and carrier tracking. An overview of GPS interference mitigation 
techniques including spatial and temporal adaptive filtering methods can be found in [29] . 

GPS ARRAY SIGNAL SIMULATOR 

In order to study various beamforming algorithms, we developed a MATLAB-based GPS 
simulator that generates signals coming from an antenna array and allows for testing of many 
different orbital scenarios and algorithms. The simulator obtains position and velocity of the GPS 
satellites from stored ephemerides and uses them to calculate the LI Coarse /Acquisition (C/A) 
GPS signals as received by each element of an antenna array. The simulator can also introduce 
near-far and multipath interference at the discretion of the user. This section provides some detail 
on the signals generated by the simulator. 

We represent the GPS LI signal at the transmitter at GPS time t as (assuming zero initial phase) 

sr(t) = V2Pb(t) coscvot, (1) 

where P is the power of the transmitted signal, ay = 2/t/q, where fo is the GPS LI 
frequency of 1.57542 GHz, and b(t) = c(t)d(t) is the real baseband signal comprised of the 
binary {+1, — l}-valued GPS C/A code signal c(t) and data message d(t). The code c(f) is a 
lms-periodic length-1023 Pseudo-Random Number (PRN) sequence designed to have strong 
auto and cross-correlation properties, which are critical for precision ranging and Code Division 
Multiple Access (CDMA). The data d(t) is a 50Hz binary data modulation whose transitions are 
synchronous with the first transition of the 1ms code [30]. 

The signal is modeled as being transmitted from the GPS satellite to a receiving antenna 
through a free-space channel model that currently ignores the effect of the ionosphere and 
atmosphere. The channel model attenuates the signal to power level C, delays it by r (t) and adds 
bandpass White Gaussian Noise (WGN) no(t) with constant spectral density No/2 over the signal 
bandwidth. 2 Here r (t) = r(t)/c where r(t) is the distance from transmitter at transmit time 
t — r (f) to receiver at receive time t, and c is the speed of light in vacuum. The signal received at 
the antenna at time t is 


s R (t) = \JlCb(t — t) cos a>o(t — t) +n 0 (t). (2) 

The model for the receiver RF front-end filters, amplifies, and down-converts the signal and noise. 
The filters pass the signal, but limit the noise to bandwidth B. A quadrature down-conversion 
process mixes the signal with the Local Oscillator (LO) generated tone e~ ,u ’ Lot to bring the signal 
and noise to a lower center Intermediate Frequency (IF) oyp = to o — a?LO, after which, low-pass 
filters remove the resulting high frequency ay + ay o components. Finally, Automatic Gain Control 
(AGC) circuits rescale the signal plus noise to an appropriate level for the Analog-to-Digital 
Converter (ADC) to quantize the data with full precision. The specific value of the rescaling factor 

2 The noise has the bandpass representation 

no(t) = n,(f) cos coot — n Q (t) sinoyf 

with baseband in-phase and quadrature components n lr n Q and n(t) = «i(t) + jn Q (t) is the complex envelope of no, 
defined by h(t) = («o(0 + /n o(t))e~> u ’° t and Ho(f) is the Hilbert transform of «o(0 [31]. In reality, this noise is primarily 
generated in the first amplifier of the receiver which also greatly amplifies the signal, effectively fixing the C/ No for all 
further processing. It is standard modeling practice to attribute this noise to the channel. 
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is inconsequential for the model, so we take it to be ^= , which conveniently normalizes the noise 

power. Defining the notation 7 = C/Nq and a = \/NoB, the resulting model at the input to the 
ADCs is 

SffO) = s R {t)2cr^ 1 e~ jcVL ° t 

= \J (2y/B)b(t — r) cosu>o(t — T)e~i c ° wt + 2cr^ 1 no(t)e~i a ’ L ° t 
= \J ( 27 /B)b(t — r)e ; ( a;iFf-a ’ oT ) + cr“ 1 u(f)e- /a?IFt 

= s(t)+n(t). (3) 

where n(t) is the complex envelope of Ho(f) (see Footnote 2 or [31]) and n(f) is then complex 
bandpass WGN noise centered at ccif with independent unit variance real and imaginary 
components. 

The signal plus noise is then sampled to get x k = x(t k ) at the regular sampling times t k = kT s + 
to, for k a positive integer, and to an arbitrary starting time. In our simulation, we assume that 
T s = 1 / B. The advantage of this assumption is that the simulated noise is then uncorrelated from 
sample to sample. Finally, we arrive at the IF sampled data model 

x k = Q[s k + n k ], 

s k = y/ (2T s7 k )b(t k - ^yW-wo**). ( 4 ) 

where Q[-] is the quantization function. 

The GPS signal simulator generates files of complex x k samples for a variety of configurable input 
parameters including the bits of precision for the quantizer Q[-], cvjp, T s , 7 (f), and r(f). The model 
(4) provides the correct IF signal structure for any GPS receiver with matching configuration 
values of cu\ R , T s , etc., to acquire and track the sampled data. 

Higher level GPS processing, e.g., solving for position, velocity and time, requires multiple 
signals from different transmitters. Each signal must have d(t ) and r(f) consistently modeled in 
accordance with GPS-IS-200 [30] to account for transmitter and receiver dynamics. The most direct 
way to model multiple transmitters consistently is to take stored GPS ephemerides and almanac 
data from Internet archives and compute transmitter states according to the algorithms outlined 
in the GPS-IS-200. 

The simulation must simultaneously reconstruct the d(t) using the ephemeris and almanac data 
along with simulation time information. Our receiver also models the receiver dynamics based on 
simple earth-orbital or random-walk dynamics. Once the transmitter /receiver (tx/rx) dynamics 
are established, the signal parameters can be calculated. Since the dynamics are very smooth the 
tx/ rx states and signal parameters r and 7 can be calculated at a relatively low rate (1 ms in our 
simulator) with simple interpolation used to determine r (t k ) and 7 (t k ) at every sampling time. 

To model the signals from a constellation of M GPS transmitters received by an N-element array, 
we can use the single tx/ rx model (4) separately M x N times. We annotate s, r and 7 with two 
indicies, i and j, to represent the path from transmitter i to array element j. The combined signal 
received by array element j can then be modeled as 

Xj = S j + Hj, 

M M . 

s i = E s a = E (5) 

i= 1 i= 1 
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where index i is summed over since the transmitted signals are received additively, and rij is 
an independent additive noise term for each array element. We have omitted the quantization 
operation and the sampling superscript k here to simplify notation. 

Two approximations can be used to simplify the model (5). The first is based on the assumption 
that the array is much smaller than the transmit-to-receive path lengths, i.e., |p ; | <C r,/, for all 
i,j, where p ; is the position offset of the j th array element relative to reference element 0. This 
assumption allows the approximation Xjj ~ t,o + k - p ( 7ay, where k, is the propagation vector 
from transmitter i to the array. Second, under the assumption that that modulation envelope 
changes negligibly in the time it takes the signal to propagate across the array, sometimes called 
the narrowband condition, we have bj ( t — r,-y ) « bj(t — r, 0 j . Together, these approximations, both 
valid in our case, give the signal at the j th element as 

M , 

Sj ~ E V (2T s 'Yi j )b i (t - Xi0 )e^-^- v i 
i= 1 V 
M 

= E s iohijr with hjj = e~i k ‘ p >. (6) 

i= 1 

Under these approximations, the signal received from the /th transmitter by each array element is 
related to that received by one reference element by a simple phase shift hjj. We refer to these phase 
shifts as the geometric phase shifts since they are defined only by the geometry of the transmitters 
and array elements. In the traditional beamformer described below, the complex conjugates of 
these phase shifts are used as weights in combining the array signals. The effect is that the signals 
are added coherently, giving maximal gain in direction of each transmitter. See Figure 3. 



Figure 3: Array geometry for a two-element array. 

In summary, the software array signal simulator, which was used for testing several beamforming 
methods, works in the following three step process: 

1. Simulate transmitter and receiver dynamics (including attitude dynamics for the latter) at 
low rate time steps (nominally 1 ms). Determine signal parameters Tjj(t), 7 jj(t), as a function 
of the tx/ rx dynamics or otherwise, as directed by input options. 

2. Interpolate signal parameters to sampling times. 
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3. Compute (5) (with approximations used in (6)) and store the result to a file. 

To simplify implementation and analysis, the following assumptions were used in the simulator. 
Transmitter and receiver clock errors, ionospheric delay, and dispersive atmospheric effects were 
not modeled. Also, to allow for easy modeling of interference and multipath, carrier-to-noise 
density ratio levels were defined by user-specified profiles rather than based on a tx/ rx link model. 

Modeling Near-Far Interference and Multipath 

One goal of this work is to demonstrate the utility of the array in mitigating two important 
interference effects: near-far interference (where the receiver attempts to process weak signals in 
the presence of much stronger GPS signals) and multipath (where receiver collocated physical 
structures provide reflection paths which cause a delayed version of the signals to appear at 
the antenna). Multipath can be one of the dominant sources of error for a GPS receiver as it can 
compromise the ability to make precise range measurements. 

Since the simulator accepts user specified C/No profiles, near-far interference can be generated 
by including weak and strong signals simultaneously. Differences higher than around 24 dB are 
known to exceed the natural cross-correlation protection of the C/A code [32], Multipath signals 
can be modeled in the simulator by replicating the transmitter signals seen by the receiver with 
additional delays, from different directions of arrival, and adding them. 

ARRAY SIGNAL PROCESSING ALGORITHMS 

The elements of an antenna array sample incident signals in the space domain like 
analog-to-digital converters sample signals in the time domain. Likewise, array signal processing 
techniques allow for the synthesis of spatially selective filters (in k space) that condition and 
extract information while reducing noise and interference, which is analogous to the synthesis 
of frequency selective filters (in to space) in Discrete-Time Signal Processing (DTSP) that do the 
same. While there are many different ways to combine the array data [33], there exists a simple, 
yet fundamental, class of algorithms that are mathematically similar to the ubiquitous N-tap Finite 
Impulse Response (FIR) filter of DTSP. These algorithms multiply the input signals Xj by complex 
weight iv j and sum over the N elements, 

N 

z = Y J wJ x j = w H x, (7) 

M 

where the asterisk * denotes complex conjugate and superscript H denotes Hermitian (complex 
conjugate transpose). Note that the array output and weight vector have been represented by 
column vectors and the index i has been removed to simplify notation. We restrict attention to this 
class of methods in our investigation. 

The beam pattern corresponding to a given weight vector w, 

N 

B( k) = X>*^ k P/, (8) 

7=1 

summarizes the response of the array to incident plane waves (exactly) and narrowband signals 
f(t ) (approximately): the response to x(t) = f(t)e^ wt ~ k ' p) is B(k)f(t)ei wt . By control of the array 
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weights, the array's beam pattern can be optimized to achieve certain desirable characteristics 
including: directivity, signal-to-noise ratio gain, and interference mitigation. This optimization 
can be done with statically fixed weight designs or with dynamically adaptive algorithms where 
the performance index evolves with the signal and interference environment over time. If the 
optimization can converge faster than the objective changes, then performance goals can likely be 
met. 

Digital vs. Analog Beamforming 

One should note that the input x in (7) could equally well be a continuous-time analog signal, or a 
discrete-time digital signal. While analog and digital beamforming can achieve the same goals in 
principle, the digital approach has the advantage that the associated front-end captures all relevant 
information available to the array. Thus, the flexibility of modern re-programmable digital logic 
and high speed computers can be used to implement and evaluate different array-processing 
algorithms. Analog beamformers have reduced computational requirements, however, and may 
be the only option for wideband signals, where digital sampling rates required for the digital 
approach may become prohibitively high. 

Traditional Beamformer 

The traditional beamformer uses the weights 

IV j = ^-e _ t k oP; (9) 

1 N 

for element j, where ko is the wave vector of the signal of interest and p j is the position of antenna 
j. Since these weights are precisely the geometric phase shifts, defined in (6), it is clear that (under 
the narrowband assumption) the different signals from the array will be combined coherently. 
Also from (8), B ( ko ) = 1. Since the signal then adds coherently, and the noise incoherently, the 
Signal-to-Noise Ratio (SNR) gain of the signal received along direction ko is N: this is the best that 
can be done in terms of pure SNR gain [33] . 

Phase compensation in the traditional beamformer is the simplest of the tested beamforming 
methods. This algorithm is not adaptive and it requires prior knowledge of the geometry of the 
array pj and the direction of the signal of interest ko- While it can steer the main beam of the array 
pattern to that direction and has reduced gain elsewhere, which can help against interference, it 
does not actively control its sidelobe or null structure. 

Optimization Based Methods 

Minimum Mean Square Error (MMSE) beamformers minimize the mean square error of the 
combined inputs from the antenna array with respect to a reference signal r. The popular Least 
Mean Squares (LMS) algorithm was implemented in the NavSDR to solve the MMSE problem 

minENr — w H x| 2 ]. (10) 

W 

In our investigation, the reference signal r was set to be the locally generated GPS carrier spread 
by the unique C/A code for each GPS satellite. The LMS algorithm creates a beam pattern that 
points toward this satellite, without requiring prior knowledge of the direction of arrival of the 
transmitted signal. 


9 



Minimum Output Energy (MOE) beamformers minimize the total output energy except in the 
direction of the transmitter. One example which was tested in simulation is the Minimum Variance 
Distortionless Response (MVDR) beamformer. The algorithm selects the weight vector w such 
that the variance of the combined signals from the antenna array is minimized, subject to the 
constraint that the signal power is fixed in the desired look direction where h is the vector of 
elements hj = Formally, it attempts to solve 

minE[|w H x| 2 ], w H h = 1, (11) 

W 

(Note, referring to (8), that w H h = B (ko).) Under this constraint, the algorithm requires knowledge 
of the direction of the signal of interest, but not of the direction of interference. 

The MVDR and LMS algorithms were implemented in the NavSDR according to the method 
described in [25]. For both methods, the step size parameter of the update equations were chosen 
by stability criteria shown in [18] and also to find the ideal balance between settling time and mean 
error. In future work, these parameters will need to be adjusted to remove their dependence on 
input signal power. 

Comparison of Beamforming Algorithms 

Traditional beamforming is the simplest of the three methods studied and has the lowest 
computational cost. However, to implement this method, and any other that explicitly points 
the array in a certain direction (including MVDR), requires knowledge of the array attitude and 
calibration data to correct for inter-element phase differences due to hardware. Table 1 summarizes 
many of the trade-offs between the three beamforming methods studied. As an example. Figure 
4 illustrates the simulated response of the LMS and MVDR methods in the presence of multipath 
interference. With the MVDR method, the main beam is centered on the signal of interest while 
the side beams are steered away from the interfering signals as much as possible for the geometry 
of this simulation. With the LMS method, in this simulation, the side beams are steered toward 
the multipath signals. In theory, there are some conditions under which LMS will mitigate the 
multipath signals, including when the SNR is lower than the SNR of the desired signal and when 
the update equation is properly tuned to the dynamics of the signal environment. 

The three beamforming methods studied were implemented in a MATLAB tracking program for 
initial verification, and later in the NavSDR. It was observed that an important modification to the 


Table 1: This table summarizes the advantages and disadvantages of implementing, in a GPS receiver, the 
three beamforming methods presented in this paper. The best method to use (highlighted in gray) varies 
depending on knowledge of the direction of the signals of interest and availability of processing power. 


Beamforming Method 

Traditional 

LMS 

MVDR 

Calculation 

direct 

iterative 

iterative 

Computation cost 

low 

high 

high 

Requires known direction of signal of interest 

yes 

no 

yes 

Requires calibration 

yes 

no 

yes 

Reduces near-far interference 

may 

yes 

yes 

Reduces interference from multipath signals 

may 

may 

yes 
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Figure 4: Simulated hemispherical beam pattern, as seen from above, of an antenna array with 16 elements, 
in the presence of multipath signals. 


LMS method was necessary due to the GPS signal structure. The LMS method requires a known 
reference signal, but since the GPS binary data is unpredictable, the algorithm was modified 
to become insensitive to the bit transitions in the data. Without this modification, the method 
becomes unstable. 

DESIGN OF ARRAY AND ARRAY FRONT-END 

An antenna array was designed and fabricated using inexpensive GPS LI patch antennas and 
characterized in one of Goddard's anechoic chambers. The array consists of 16 elements with 
equal horizontal and vertical spacing between antennas, in a 4x4 grid. The spacing is half of the 
GPS LI wavelength (A/| / 2), about 9.5 cm, which is optimal for a regular array with a beam pattern 
having a single, narrow-as-possible, main beam. The main beam width is inversely proportional 
to the size of the array and wider spacing than An / 2 leads to multiple main beams or "grating 
lobes" (analogous to aliasing in DTSP for sampling rates lower than the Nyquist rate). 

An array RF front-end was also designed and built for laboratory use, to coherently down-convert 
and sample each of the signals from the individual array elements. The design is partitioned 
into a motherboard and 16 daughterboards (each one producing digital samples from one of 
the antennas). The samples are filtered and packed in an FPGA on the motherboard, which 
also handles frequency synthesis and hosts the power circuits. Figure 5 shows a diagram of the 
front-end, and Figure 6 shows an image of the complete motherboard with daughterboards after 
they were built and assembled. 

ARRAY RECEIVER TEST RESULTS 

The three beamforming methods were implemented in the NavSDR and tested using both the 
MATLAB simulated and live-sky signals. In all the live-sky tests our 16-element antenna array 
and custom front-end were used. Tests of multipath and near-far interference were conducted 
using the MATLAB simulator since such scenarios are impossible to control with live-sky signals. 
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Figure 5: Block diagram of array front-end design. 



Figure 6: Array front-end motherboard with 16 Figure 7: Array front-end single-channel 

daughterboards (partially visible) attached downconversion and digitization daughterboard 


Simulated Signals 

Single Antenna and Array (Without Interference) 

The first test simulated GPS signals received from five satellites, at a fixed C/No = 45dB-Hz as 
seen by a single antenna. This level is typical for a receiver on the Earth's surface. Figure 10a shows 
the array-receiver estimates of the C/No for the five satellites, first using one antenna and then, 
after 20 seconds, combining all 16 antennas using the LMS beamforming method. The observed 
increase in C/No was approximately 12 dB, expected for 16 antennas (as 101og 10 16 = 12 dB). 
The results using the other two methods showed similar C/No increases. We also studied the 
effect of the array on the Position, Velocity, Time (PVT) solution quality. The standard deviation 
of the calculated position was improved from about 14.5 m before beamforming to about 4.5 m 
afterwards. This error is somewhat higher than expected for the corresponding signal levels and 
geometry (8.8m and 2.2m, respectively); the reason for this has not yet been identified. Regardless, 
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the fractional decrease in the error (69%) is close to the expected improvement (75%). 

Near-Far Interference Test 

To simulate near-far interference, four GPS transmitters were set at a signal level of 35 dB-Hz, 
and a fifth at 65 dB-Hz. The results in Figure 8 show the C/No for several GPS transmitters in the 
presence of an interferer. In each scenario, the interference starts 20 seconds into the simulation. 
Figure 8a shows the signals when they are tracked using one antenna. Once the interference 
begins, the signals are severely attenuated and several cross-correlations are tracked, causing the 
PVT solution to be corrupted. In Figure 8b, all the signals from the array are combined using 
the traditional beamforming method. The gain in the direction of the signals of interest provide 
some protection against the interferer. The LMS method (Figure 8c) fails in this case. This is not 
unexpected because the interference exceeds the 24 dB cross-correlation protection level. The 
interferer and local reference have a strong enough correlation so that the algorithm converges 
on the interference and not the desired signal. The MVDR method (Figure 8d) improves the gain 
in the direction of the signal of interest (over the traditional beamformer), as well as reducing 
interference. 

Multipath Test 

In this scenario, signals from five GPS transmitters were simulated. For each transmitter, the 
direct signal and one multipath signal were created, the multipath signal being a replica of the 
corresponding direct signal with a different delay (20, 40, 60, 80 or 100 meters) and direction of 
arrival. All signals were simulated with the same C/No of 45 dB-Hz, representing an extreme 
case of multipath interference. Performance was measured by the effect on the point solution 
position. As shown in Figure 9, with one antenna, the multipath interference causes a large 
bias and increased noise level in the calculated position. This result exemplifies the effect of 
multipath interference in GPS performance: the multipath signals distort the code correlation 
function, causing bias in the measurements, while also causing destructive interference of the 
carrier, reducing the C/No ■ With any beamforming method, the C/No is increased, which results in 
a better PVT solution due to less noise, but not necessarily a lower overall error as the position may 
still be biased. In this simulation, only the methods based on direction of arrival (the traditional 
beamformer and MVDR) are able to attenuate the effect of the multipath signals significantly and 
reduce the overall error. 

Live-Sky Signal Tests 

Array Calibration 

The beamforming methods were also tested with live GPS signals, using our low-cost antenna 
array. While the simulated data assumes that the phase differences between the antennas are only 
due to their position and the direction of arrival of the signals, the RF front-end has unequal path 
lengths, causing additional phase differences between each of the signals from the antennas. These 
phases have to be calibrated out for the beamforming methods that require a known direction of 
arrival (traditional beamformer and MVDR). The calibration was done by tracking the signal from 
the same satellite with all antennas and comparing the phase of the local carrier calculated by the 
tracking loop of the receiver software. The procedure was repeated for several satellites, and the 
phase difference due to position of the antennas and direction of arrival was subtracted from 
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(a) Without beamforming 



(c) LMS 



time (s) 

(b) Traditional beamforming 



Figure 8: C/Nq before and after strong near-far interference begins (20 seconds into the simulation), with 
and without beamforming. Each trace represents the signal from a different satellite. The interference 
causes a large decrease in C/Nq with one antenna (subfigure 8a) and the acquisition and tracking of 
cross-correlated signals. The 65 dB-Hz interferer is not tracked because the signal power is too high for 
the NavSDR (certain calculations overflow and fail). Traditional beamforming (8b) and MVDR (8d) reduce 
the effect of the interference. In each case, there is only a small decrease in C/ No and cross-correlated signals 
are avoided. LMS (8c) is not able to maintain track of any signal in this case. 


the total phase difference from the local carrier. The remainder was approximately constant, and 
assumed to be the bias. 

Single Antenna and Array (Without Interference) 

The receiver test results in Figure 10b show C/No estimates for five satellites initially using only 
one antenna, and after 60 s, combining the signals from all of the antennas using the LMS method. 
The increase in C/No upon turning on the array was about 11 dB on average, close to the expected 
12 dB of gain for 16 antennas. Unlike in the simulations, the increase in C/No varied depending 
on the satellite, probably due to the varying sensitivities of the antennas for different directions of 
arrival. Also, the initial C/Nq estimates are somewhat lower than expected for live-sky signals. 
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one antenna 
traditional 
LMS 
MVDR 



Figure 9: Error in position calculation in the presence of multipath interference. Multipath signals start 
at about 17 seconds. LMS reduces the variance in the calculated position, but not the bias. Traditional 
beamforming and MVDR reduce the overall error caused by multipath signals. 


This loss is attributed to issues with the array front-end which are still under investigation. 
As expected, the error in the position calculated by the receiver decreased with the increase in 
C/No. The standard deviation of the calculated position was improved from about 14 m before 
beamforming to about 5 m afterwards. These levels are consistent with our earlier investigation 
using simulated signals, presented above. 



(a) Simulated signals 



Figure 10: C/No of GPS signals before and after beamforming. The increase is close to the ideal value (12 
dB for 16 antennas). 


CONCLUSIONS AND FUTURE WORK 

In this paper, we presented preliminary results from the Navigator GPS Digital Beam Former 
(NavDBF) technology development project. The paper described the design and development of 
a 16-element antenna array and front-end electronics, the development of a MATLAB-based GPS 
signal simulator, the evaluation and implementation of one fixed and two adaptive algorithms 
for Digital Beam Forming (DBF) in both a simple off-line MATLAB GPS receiver simulation and 
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our NavSDR real-time software GPS receiver development platform, and finally the assembly of 
a complete GPS receiver with beamforming capability Promising initial results were presented, 
using both simulated and live-sky GPS signals, which largely demonstrated the expected gains in 
sensitivity and interference mitigation for both near-far and multipath interference. 

Future work will focus on pushing the NavDBF along the standard three-stage technology 
development process, which was also introduced in this paper. The next stage in the process 
involves porting the algorithms from the NavSDR to the NavDev FPGA platform. We are also 
developing a flight-qualified, miniaturized, RF front-end. A GSFC IRAD proposal was selected to 
begin this work, which will leverage a related NASA Phase II Small Business Innovation Research 
(SBIR) contract that is funding the development of a multichannel GPS front-end Application 
Specific Integrated Circuit (ASIC). 

Future Applications for NavDBF 

The NavDBF has potential application to missions whose attitude dynamics do not allow for a 
permanent zenith (for LEO) or nadir (for GEO) pointed GPS antenna, and who may, therefore, 
have some difficulty maintaining good GPS coverage with fixed body-mounted antennas. It also 
has application to missions that need high precision in PVT measurements in a poor multipath 
environment such as on the International Space Station (ISS), specifically the multipath mitigation 
benefit. 

The MMS mission, in its highly-eccentric, high-altitude orbits, and use of four spacecraft in 
formation, represented an ideal target application for the Navigator GPS receiver. In addition 
to the weak signal environment, physical constraints on the spacecraft design forced the use of 
an array of four GPS antennas placed around the perimeter of the spacecraft. As the four MMS 
spacecraft rotate (at 3 revolutions per minute), the receiver maintains continuous GPS tracking 
by switching from one antenna to the next in sequence. This switching mechanism works well 
enough to meet mission requirements but has performance costs relative to an omnidirectional 
antenna. A beamforming approach, which would combine the array of signals, would have more 
gain and allow for smoother tracking of the signals during spacecraft rotation. This approach 
could be applied to future missions with similar challenges. 

An interesting example application of NavDBF is a hypothetical Geosynchronous Earth Orbit 
(GEO) spacecraft-servicing mission that would need to maneuver around a GEO target spacecraft 
and may not be able to support a nadir-pointed or gimbaled high-gain GPS antenna. Such a 
mission would benefit greatly from an array of GPS antennas that could be optimally combined. 
For such a mission, the potential application of NavDBF goes beyond just providing effective 
GPS navigation. The NavDBF could provide a sensor capability for Automated Rendezvous and 
Docking (AR&D) operations by tracking GPS signals reflected from the target spacecraft and 
measuring their angle of arrival and range to target. The Navigator team has investigated this type 
of application during HST-SM4 and is currently actively involved in continuing this research. 
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Introduction 



• Global Positioning System (GPS) 

- Satellite navigation system developed by the United States military, 
became fully functional in 1995 

- Signal levels and availability are adequate for users on or near Earth 

• Navigator 

- GPS receiver developed at NASA for use in high orbits 

- Increase in sensitivity achieved with fast acquisition 

- Ten times more sensitive than traditional GPS receivers (25 vs. 35 
dBHz) 

- Currently included in Global Precipitation Measurement (GPM) and 
Magnetospheric Multiscale (MMS) missions 


Next-generation Navigator goals 

- Further increase signal-to-noise ratio (SNR) by 10 dB or more (15 dBHz) 

- Reduction of interference in near-far problem and multipath signals 
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Traditional GPS receiver: 
sensitive to 35 dBHz 
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Current-stage Navigator: 


sensitive to 25 dBHz 
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Next-generation Navigator: 

sensitive to 15 dBHz 
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Navigator: interference environment 



Near-far problem: strong signal from a 
nearby transmitter prevents acquisition of 
weaker signals from transmitters farther 
away 

- Occurs when receiver is close to a GPS 
satellite 



strong signal 
receiver 

weak signal 


transmitter 
(satellite) 
Near-far problem 


Multipath interference: combination of multiple rays degrades signal 
level and causes error in GPS measurements due to additional 
delays 

- Occurs near reflective surfaces, common on Earth and in spacecraft 

Ins, it is the 
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(satellite) 


direct 
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receiver 
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Antenna arrays and beamforming 




• The use of an antenna array with beamforming methods is an 
effective way to achieve both an increase in SNR and reduction of 
interference 

- Beamforming is the formation of beam patterns through the selection of 
coefficients or excitations to each antenna to combine the signals 

- Ideally, the main beam is steered toward the signal of interest, 
amplifying it, and nulls are steered toward interference, reducing them 
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Navigator with antenna array: components 



Antenna array: 

• signal reception 

• amplification 


Front-end: 

• frequency 
conversion 

• filter 

• sampling 



Navigator GPS receiver with antenna 

array 

Software: 

• simulations 

• beamforming 

• acquisition and tracking of GPS signals 

• calculation of position, velocity and time 
(PVT) 
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GPS signal simulator 



• Simulator needed to study beamforming algorithms and test 
different signal levels and interference scenarios 

- Calculates dynamics of satellites and receiver, in low time steps 

- Generates GPS signals consistent with dynamics, interpolating to 
sampling times 

- Computes signals from each satellite as received by each antenna in an 
arbitrary array, approximating delays between antennas as phase shifts 

- Stores simulated signals in a file, to be read by the Navigator GPS 
receiver software 

- Supports near-far and multipath interference 


dynamics 



GPS signals antenna array 





interference 

(optional) 



GPS signal simulator sequence 


GPS signal file 
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1101000111 

0100101011 

0101010101 

0100101010 

1010101010 

1001010101 

1011010101 
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Beamforming algorithms 





Beamforming method 

Traditional 

Least mean squares 
(LMS) 

Minimum variance 
distortionless 
response (MVDR) 

Criteria 

maximize gain at 
desired direction 

/k r, 

Wj = Q J 1 

minimize difference 
from reference signal 

. F r q 

mmE r-^j^Xi 

w 

minimize power except 
at des red direction 

• 2 1 

mmE h=l W i X i 

w J 

-/kr, i 

E,=i«i e 1=1 

Calculation 

direct 

iterative 

iterative 

Computation cost 

low 

high 

high 

Requires known direction 
of signal of interest 

yes 

no 

yes 

Requires calibration 

yes 

no 

yes 

Reduces near-far 
interference 

may 

partially 

yes 

Reduces multipath 
interference 

may 

no 

yes 


The best method to use depends on knowledge of directions of 
signals and availability of processing power 
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Results: simulated and live-sky signals 

(without interference) 
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five independent beam 
patterns formed, one for 
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Results: simulation of near-far interference 

interferer at 65 dBHz, after 20 s 
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Results: simulation of multipath interference 

I SPACE FLKJMT CENTER 

• One direct signal and three multipath signals, for one GPS satellite 

• All signals have equivalent power and differ only in phase 

• signal of interest 

Beam pattern, LMS Beam pattern, MVDR ♦ interfering signal 



Beam pattern of antenna array under multipath interference 
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Results: simulation of multipath interference 



• One direct signal and one multipath signal, for each of five GPS 
satellites 

• Each multipath signal has a different delay and direction of arrival 

• All signals have equivalent power (45 dBHz) 
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Conclusion and future work 



• A digitally-steered antenna array was simulated, designed and built 
for the Navigator GPS receiver 

• Preliminary results from both simulated and live-sky signals are 
promising, and demonstrate the expected gains in sensitivity and 
mitigation of interference 


• Future work: mature this technology and prepare for flight 

- More thorough analysis of interference scenarios and stability of 
beamforming algorithms 

- Port beamforming capability to Navigator platform in a field- 
programmable gate array (FPGA) and miniaturized front-end 

- Possible applications in satellite-servicing missions, such as Restore, 
and in experiments at the International Space Station (ISS) which 
require precise position or timing, such as the Neutron Star Interior 
Composition Explorer (NICER) 
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